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I. INTRODUCTION 
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The probability representation, in which cosmological quantum states are described by a stan- 
dard positive probability distribution, is constructed for minisuperspace models selected by Noether 
symmetries. In such a case, the tomographic probability distribution provides the classical evolution 
for the models and can be considered an approach to select "observable" universes. Some specific 
examples, derived from Extended Theories of Gravity, are worked out. We discuss also how to 
connect tomograms, symmetries and cosmological parameters. 
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The main issue of Quantum Cosmology is to find out methods, derived from some quantum gravity approach, to 
^vq ■ achieve the so called Wave Function of the Universe (or the corresponding density matrix), a quantity related to the 
initial conditions from which dynamical systems (representing " classical universes" ) can emerge 0, • 

Different representations of the density matrix such as the Wigner quasi-distribution were applied to study cos- 
mological models in the minisuperspace (all- Recently, a tomographic probability representation of quantum states 
was suggested: the quantum state can be associated to a standard positive probability distribution (tomogram). 
The tomographic probability of quantum states was introduced in Quantum Cosmology in 5]. In this framework, 
__i | cosmological models described by the solutions of the Wheeler-de Witt equation can be mapped onto solutions of the 
Fokker-Planck-likc equations for standard probability distributions. The tomographic probabilities are connected to 
the Wigner function of the Universe by the Radon integral transform 0] and by the wave function of the Universe 
with the wave related to the fractional Fourier transform [H, [f| considered in quantum mechanics in @, Q • 

Some important remarks are in order to point out why the tomographic approach is useful to improve the Wigner 
picture. As we said, the Wigner function [9( can be used to describe quantum states which are "universes". It has 
£f) • been introduced in order to deal with mixed quantum states with some remarkable properties. In contrast to the 
standard wave function and density matrix representation, where complex functions are used, the Wigner function is a 
real one. It depends on two real variables, p and q, which label a point in the phase space. The marginals constructed 
' by the Wigner function W(p, q) are integrals either over q or over p. They are positively defined probability densities 
f^) . in position and in momentum respectively. These properties make the Wigner function to be very similar to the 
probability distribution density used to describe the classical states of a particle in the phase space as in classical 
statistical mechanics. The Wigner function was introduced, namely, to find a rigorous description of quantum states 
in terms of objects which have a classical probability distribution. However, such a function cannot play this role since 
it can take negative values and this feature does not allow to interpret the function as a fair probability distribution 
density. The necessity to find out such a kind of probability density, which rigorously describes the quantum state, 
is known as the Pauli problem 10]. Pauli supposed that the probability distribution in momentum, together with 
the one in position, determines completely the quantum state. But it was shown that different quantum states can 
exist having the same probability distribution density, so the Pauli problem was proved to have no solution, if one has 
only these two probability distributions (namely the one in momentum and the one in position). Nevertheless, this 
problem got a complete solution in the tomographic probability approach to quantum state description. In short, one 
needs a family of marginal probability densities in the phase space obtained from rotations in phase space of the Pauli 
distribution pair: such rotations can be parameterized by a given angle 8. This family, which is nothing else but the 
Radon integral transform [ll[ of the Wigner function, is sufficient to completely describe the quantum state, i.e. to 
find the density matrix of Wigner function if one knows the tomographic probability distribution. The tomographic 
probability contains the same information on the quantum state as the Wigner function of density matrix. On the 
other hand, like the Wigner function, being completely equivalent to density matrix in position and momentum 
representations, the tomographic probability density has its own merit: it is a fair "positive" probability distribution. 

It is well known that the Wigner function has been intensively used in Quantum Cosmology [12j |. Besides, some 
applications of tomographic probability, based on its property to be the fair positive measurable probability distribu- 
tion, have been considered for cosmological problems [a, |q| . For example, the notion of tomographic entropy, using 
the Shannon construction of entropy and information [13 |. provides extra characteristics of cosmological quantum 
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states (g, ||| . The Shannon entropy cannot be used for the Wigner function since such a function is not a probability 
distribution. 

Furthermore, using the properties of tomographic probability, related to its own features, it is possible to address 
several problems of Quantum Cosmology whose solutions cannot be achieved by the Wigner function. In particular, 
it has been shown that symmetries can play an important role to select viable cosmological models [lj, EH Ell- 
in the framework of the so-called Noether Symmetry Approach (l6l . ITtI EH it is possible, in principle, to find out 
cyclic variables related to conserved quantities and then to reduce cosmological dynamics. Besides, the existence 
of symmetries fixes the forms of couplings and self-interaction potential giving the relation between them in the 
interaction Lagrangians. The existence of Noether symmetries for minisuperspace cosmological models can be viewed 
as a sort of selection rule to recover classical behaviors in cosmic evolution In fact, it was shown in [l[ that the 
form of the wave function with specific picks provides a classical regime in the evolution of the Universe. The picks 
(as oscillations) are characteristic properties of such regime and are selected in cosmological models where Noether 
symmetries are present (l9j . 

In the framework of minisuperspace approach to Quantum Cosmology, the aim of this paper is to study cosmological 
models in order to seek for conditions to select classical (observable) behaviors by tomographic probability represen- 
tation. We find out that, according to 19], if tomograms are related to extra symmetries (Noether symmetries), 
they allow oscillatory behaviours of the wave function and then give rise to classical solutions of dynamics. In the 
framework of the Hartle criterion, such solutions are observable universes. The approach is worked out for general 
classes of Extended Theories of Gravity (see for example [20|, l2l| which have acquired a huge interest in recent years 
as possible schemes capable of explaining cosmological dynamics at all epochs, starting from inflation, through matter 
dominated era, up to the today observed accelerated behavior [22|, HH • 

The paper is organized as follows. In Sec. 2, the tomographic approach and its relation to the Wave Function of 
the Universe, considering the Wigner function and density is reviewed. Sec. 3 is devoted to the semiclassical limit of 
Quantum Cosmology and the Hartle criterion to select observable universes. In Sec. 4, the stationary phase method is 
discussed in view of the WKB approximation for the Wave Function of the Universe. Sec. 5 is a short summary of the 
minisuperspace cosmological models, coming from Extended Theories of Gravity, which we are going to analyze in 
the tomographic representation. Exact solutions (wave function of the universe) are obtained if Noether symmetries 
exist. Such solutions show oscillatory behaviors and then the possibility to implement the Hartle criterion. In Sec. 6, 
we give the related tomograms in stationary phase approximation for the previous selected solutions. Conclusions and 
a discussion of possible relations to the observed universe the classical regime evolution related to initial are presented 
in Sec. 7. 



II. TOMOGRAMS OF COSMOLOGICAL QUANTUM STATES 



We consider here the tomographic map of the wave function of the universe (or its density matrix) in the framework 
of the minisuperspace approach. Let us take into account first the minisuperspace model in which the pure space of the 
universe is described by a vector |\& > | in the Hilbert space of quantum states. The wave function ^(x) =< x\^f >, 
x G 3? can be mapped onto a fair probability distribution W(X, /i, v) called tomogram 
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This probability distribution is positive and normalized for all the parameters fj,, v G 3? 

J W{X,fi,v)dX = 1, 

if the wave function is normalized, i.e. 



I*G/)I S 
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The tomographic map of the wave function of the universe onto the universe tomogram 



(1) 



(2) 



(3) 



(4) 



is invertible (up to a constant phase factor). One has the inverse formula determining the density matrix of the 
universe pure state in the positive representation 
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Thus one has the relation 



|*(0)| 2 = — / WfafrO^dydn 

In view of Eq.©, one has the connection of the wave function of the universe with its tomogram 

1 



$(x) 



2tt|*(0)| 



W(y,n,x)e^ y -" x ^dydn 



or in view of Eq. ^ 



(5) 



(6) 



(7) 



■9i(x) = y ( yV(y,fi,x)e i( - y - fiX ^dyd f j, 



2tt 



W(2/',/x',0)e^d 2 /'^' 



-1/2 



(8) 



Eq.© determines the wave function of the universe in terms of the universe tomogram (up to a constant phase factor). 
This constant phase factor of the wave function is not essential. In fact, the wave function can be written in form of 
a modulus and a phase factor 



*(z) = \^(x)\e i3 ^\ 



(9) 



In view of Eq. |8|) one can obtain both the amplitude and the phase of the wave function in terms of the tomogram. 
For the amplitude, one has 



l*W| 2 = ±f W(y,n,x)e^-^dydfi. 



(10) 



Since the amplitude squared is a real (positive) function, it is expressed in terms of the real part of the right hand 
side of Eq.([TU|). i.e. 



I*0z)| 2 = ^ J W(y,n,x) cos(y- nx/2)dydfi. 



(11) 



We point out that the tomogram W(y, fi, x) is a real non negative function. On the other hand, the phase of the wave 
function S(x) is also determined by the universe tomogram since 



S(x) = —i In 



one has 



S(x) 



-i In 



J W{y,ii,x)e i{ -y-^ x ^dyd^ 
| J W{y, n, x)e^y-^ x /^dy dfi\ 



(12) 



(13) 



In quasi classical approximation (WKB) for the wave function of the universe, which corresponds to a stationary state 
with energy E, the amplitude and the phase are connected. In this case the phase S(x) satisfies the Hamilton- Jacobi 
equation. One can consider the tomogram given by Eq.([T]) as follows 



2ttM 



\^(y)\e lS(y)+ ^ y2 '^ v dy 



(14) 



The approach outlined above is in complete agreement with the semiclassical limit of Quantum Cosmology and it 
provides a physical interpretation of the wave function of the universe, as we are going to discuss below. 
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III. THE SEMICLASSICAL LIMIT OF QUANTUM COSMOLOGY AND THE HARTLE CRITERION 

Scmiclassical limit occupies an important role in Quantum Gravity due to the lack of a self-consistent and definitive 
theory. Only thanks to such a limit, we can obtain the solution of the Wheeler-De Witt (WDW) equation, the 
wave function of the universe VP, and search for interpretative criteria which allow physically motivated previsions. 
Furthermore, cosmological observations are based on a "classical universe" since we are not able, at least considering 
today's facilities, to obtain data toward and before the Planck epoch. Taking into account also the fact that the 
main goal of Quantum Cosmology is to give self-consistent laws for initial conditions, the correct identification of a 
semiclassical limit is crucial. In standard quantum mechanics, the wave function can be developed as a power law 
series in h in the WKB approximation. Formally, the semiclassical limit corresponds to the limit H — ► 0. Analogously, 
we can write down the WDW equation as 



-m 2 P u)*[h ij (x),(t>(x)}=0 (15) 

where U(hij,(j>) is the superpotential, hij are the components of the spatial three-metric, the geometrodynamical 
variables, and cf> is a generic scalar field describing the matter content. The semiclassical limit can be achieved by 
using nip 2 , the squared Planck mass, as an expansion parameter. 
The wave function can be represented as 

Vlhijix),^)}-^? 3 (16) 

where S is the an action. A state with classical correlations has to be a superposition of states of the form (JTHJ). The 
WKB approximation is achieved by the expansion 

S = S + m p 2 S l + 0(m p 4 ) (17) 

and inserting it into the WDW equation. Then equating terms of the same order in nip, we obtain the Hamilton- Jacobi 
equation at any order in S. At the lowest order in So, we get 

VS ■ VS + U = . (18) 

In the same way, we obtain equations at higher orders which can be solved taking into account results from the 
previous orders. It can be shown that we need only VSo to recover the semi-classical limit of Quantum Cosmology 
[24| , that is the wave function of the Universe is given by 

If So is a real number, we have oscillating WKB modes and ^ is peaked on a phase-space region defined by the 
equations 



n = m|,^ (21) 

where ELy and 11^ are respectively the classical momenta conjugate to h* 3 and <fr. 

Using the Hamilton- Jacobi Eq. (fT8|) . Eqs. (|20|) and (121|) are first integrals of classical field equations defining a set of 
solutions. This fact, in principle, can solve the problem of initial conditions. In fact, for a given action So, solutions 
(120p and (|2ip involve n arbitrary parameters, while the general solution of the field equations involves {2n — 1} 
parameters. Then the wave function is peaked (or oscillates, being So real) on a subset of the general solution. In 
other words, the boundary conditions on the wave function imply initial conditions for the classical solution. Besides, 
if the wave function 'J is sufficiently peaked, oscillating about some region of the configuration space Q = , </>}, it 
is possible to find correlations among the observables that characterize this region. If it is not peaked, correlations 
among the observables are not present. 

At this point a discussion on what we mean as "classical universe" and how the concept is related to an oscillating 
\E' is due. 
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Let us take as an example a cosmological model characterized by a set of observable parameters as Hq, the Hubble 
constant, qo, the deceleration parameter, Q^i the density parameter for matter, the density parameter for the 
cosmological constant A (or dark energy) and fi^, the amount of spatial curvature 1 . All of them are related by some 
cosmological model, solution of the Einstein-Friedmann equations and can be measured with a certain accuracy, by 
observations (see e.g. [25[ and references therein for a discussion). A "good" wave function of the universe is that one 
which is "peaked" about such a solution of the cosmological equations (which can be alternatively derived from the 
momenta (pZO)) and (|2ip and then "correlates" the set of observables {Ho, qo, SIm, ^a, ^k}- It is worth noticing that 

does not predict any specific value of the set of observables, but guarantees the correlations among them. This is 
the so-called Hartle criterion [l[ . On the other hand, such a correlation allows to define a parameter able to label the 
points along the classical trajectories, solutions of and HJ). Defining a tangent vector in the configuration space 
of classical trajectories, where * is peaked, it is 

4- = 2V5 • V (22) 

(XT 

where r is the "proper time" along the classical trajectories. From these considerations, we can state that 

1. time emerges as a parameter which labels points along trajectories where the wave function is peaked; 

2. time parametrization invariance (a symmetry of classical theory) comes out as the freedom to choose such a 
parameter. For example, assuming the redshift z, in principle, every classical cosmological solution can be 
expressed in terms of z; 

3. classical time, and in general, classical spacetime are proper notions which emerge only in a superspace regions 
where the wave function is oscillating; 

4. the existence of oscillating wave function, and then of classical spacetimes, depends on the boundary conditions 
of the theory, in particular on the shape of superpotential U(hij, <$). 

At this point it is clear what is the semiclassical limit of the theory. It is the superspace region where the wave 
function oscillates with big values of the phase which indicates strong correlation among the dynamical variables of 
the form (|2"D)) and (|2"Tj) . Such variables describe the classical trajectories in the superspace and make the concept of 
classical spacetime emerge. The tangent vector to this classical paths defines the proper time. Formally, in analogy 
with standard quantum mechanics, the semiclassical limit is obtained for mp — > oo, which means low energies with 
respect to the Planck scale. Furthermore, the limiting conditions for the wave function select a particular set of 
classical universes and the "measure" defined by the ^ itself indicates a "typical" universe. In principle, this is the 
way in which Quantum Cosmology approaches the problem of initial conditions. 

The semiclassical region of superspace, defined by the oscillating structure of this wave function, is called the 
"Lorentzian" region. This can be identified as the " ensemble" of all 3-geometries embedded in a classical spacetime. 
The region outside the Lorentzian one is called "Euclidean" . Here the action is imaginary, i.e. So — il and the wave 
function is exponential 

* ~ e -"^ 2/ . (23) 

The wave function of this form is not "classical" since it corresponds to an Euclidean spacetime. 

If W is a WKB solution , / is the action of Euclidean solutions of field equations called "instantons" . Unlike the 
Lorentzian case, the wave function (|23p is not peaked on a set of instantons. It is not "classical" since it cannot predict 
classical correlations among Lorentzian momenta and their conjugate variables. With these considerations in mind, 
it is clear that any method by which conserved quantities, corresponding to conserved momenta, can be achieved is 
useful to define a semiclassical limit in Quantum Cosmology. I other words, conserved momenta like = Eo or 
= Ei allow, in principle, to select classical trajectories corresponding to observable universes. 

In general, the Hamiltonian constraint gives the WDW equation, so that if |^ > is a state of the system (i.e. the 
wave function of the universe), dynamics is given by 

H\V>=0. (24) 



1 Actually, one should consider also the age of the universe to, but, usually, it is referred as to ~ H Q 
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In [191 ], it is shown that if Noether symmetries exist, a reduction procedure of dynamics can be implemented and 
then, we get 

-idi\V > = Ei|* > , 

-id 2 \^ > = £ 2 |tf > , (25) 



which are translations along the q J axes of configuration space (minisuperspace) singled out by the corresponding 
symmetry. Eqs. (|25[) can be immediately integrated and, being Sj real constants, we obtain oscillatory behaviors for 
I* > in the directions of symmetries, i.e. 

m 

|* >= ^ sQi \x{Q l ) > , m<l<n, (26) 
j=i 

where m is the number of symmetries, I are the directions where symmetries do not exist, n is the total dimension 
of minisuperspace. Viceversa, dynamics given by ([24]) can be reduced by (|25|) if and only if it is possible to define 
constant conjugate momenta, that is oscillatory behaviors of a subset of solutions I* > exist only if Noether symmetry 
exists for dynamics. 

The m symmetries give first integrals of motion and then the possibility to select classical trajectories. In one and 
two-dimensional minisuperspaces, the existence of a Noether symmetry allows the complete solution of the problem 
and to get the full semi-classical limit of Quantum Cosmology. As we shall see below, the tomographic probability 
representation gives results strictly related to this method since quantum tomograms become "classical" as soon as 
Noether symmetries exist. This occurrence allows, in principle, to trace back the cosmic evolution from quantum to 
classically observable states for several Extended Theories of Gravity. 

It is worth noticing that the present approach has to be compared with previous results where minisuperspace 
method has been applied to similar models. For example, in [26| . it is studied the quantum-to-classical transition in 
Jordan-Brans-Dicke quantum gravity showing the relevance of the method for inflation. In [2?], [28| the scale factor 
duality and the cosmological constant problem are faced from the same minisuperspace point of view. As we are going 
to demonstrate, the tomographic approach is fully coherent with these results. 



IV. THE STATIONARY PHASE METHOD 



Let us now evaluate the integral (| 14[) using the stationary phase approximation in view to select classical state by 
tomographic representation. The method is connected with the evaluation of an integral of the form 



A{z)e l ^ z) dz 

where z is a complex variable and a series decomposition of the function 3>(z) is used near its extremum 

i 

$(z) «$(*„) + - — {z-z ) 2 . 



The point(s) zq is(are) the solutions of the equation 



dz 



2 dz 2 



(zo) = . 



(27) 



(28) 



(29) 



The function A(z) is considered slowly varying near the point zq. Then the integral (|27[) reads 



I~A(z ) 



^/-£*"(*o) 



i*(z ) 



(30) 



In the case of N points zq — > Zq (with k = 1 , . . . N) , the integral assumes the form of a sum over these points. 

Let us consider, in an explicit exact expression (|14[) . the modulus to be a slowly varying amplitude function 

and the function 



*(g) = S(q) 



M 2 Xq 
2v v 



(31) 
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The point yo is determined by the expression 



OS fi 

7T+ 1 
aq v 



One can see that using the notation 



dS_ 



X 
v 



P 



the relation (|3"2")l is equivalent to 

fiq + up = X 

One can remind that in classical mechanics the tomogram is defined as (29j 

W{X,li,v) = J f(q,p)5(X - fiq-up)dqdp. 

In quantum mechanics, the tomogram is determined via the Wigner function W(q,p) as [3C 

dq dp 



W{X, n,v)= / W(q,p) 5(X - m - up)- 



2n 



Both relations mean the validity of The Wigner function is determined by the wave function 

W(q,p) = J * (q + |) (g _ |) e -»P du . 



(32) 



(33) 



(34) 



(35) 



(36) 



(37) 



In WKB approximation, the expression for the Wigner function is given just by the stationary phase method of 
evaluating the integral (f3"T|) that is 



W(q,p) = 1 1* (g+ |)| 



JS(l+§) e -iS{q-^) e -iup du _ 



(38) 



This expression was studied and applied in the cosmological context [!, 0| • 

Since the Wigner function determines the tomograms by (|36| . one can get the quasi classical approximation for 
the tomogram inserting into ([36| the Berry's expression for the Wigner function. But as it is clear, one can get the 
approximation directly using in (|30[) expression the function <&(<?) given by (|31[) . Then the point go is defined by 
and $ (qo) reads 



Thus one has 



$ (<?o) 



W(X, M ^)«^-|v|/(g )| 2 



d 2 S(g) 



9=90 V 



d 2 S r 

w {qo)+ ~ 



l*(9o)r 



fjr(?o)f + A< 



(39) 



(40) 



(i) 



More generally, if Eq. (|32|) has more solutions 



W(X,/i,^) 



90^' i the tomogram takes the form 

2 



N 

E 



"S9 5 " 



9o 



(41) 



The relation with the semiclassical limit of quantum cosmology is straightforward. In fact, the existence of such 
stationary points means that the wave function of the universe is peaked on conserved momenta and then classical 
trajectories in the minisuperspace can be found out. 
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V. MINISUPERSPACE MODELS FROM EXTENDED THEORIES OF GRAVITY 



The previous discussion can be applied to several classes of minisuperspaces, in particular one can take in to 
account Extended Theories of Gravity as the scalar-tensor theories or higher-order theories of gravity. Such theories 
are interesting since they are directly related to the issue of recovering suitable effective actions in quantum gravity 
[3ll | . Starting from pioneering works of Sakharov [32| . the effects of vacuum polarization on the gravitational constant, 
i.e. the fact that gravitational constant can be induced by vacuum polarization, have been extensively investigated. 
All these attempts led to take into account gravitational actions extended beyond the simple Hilbert-Einstein action 
of General Relativity (GR) which is linear in the Ricci scalar R. The motivation is to investigate alternative theories 
in order to cure the shortcomings of GR, essentially due to the emergence of singularities in high-curvature regimes. 

The Brans-Dicke approach is one of this attempt which, asking for dynamically inducing the gravitational coupling 
by a scalar field, is more coherent with the Mach principle requests [33j |. 

Besides, it has been realized that corrective terms are inescapable if we want to obtain the effective action of quantum 
gravity at scales close to the Planck length (see e.g. [13]). In other words, it seems that, in order to construct a renor- 
malizable theory of gravity, we need higher-order terms of curvature invariants such as R 2 , R^R,,,,,, R^ va P 'R uvn n, ROR, 
RO k R or nonminimally coupled terms between scalar fields and geometry as ip 2 R (see [35ll36[l37j for a review). 

Stelle [38[ constructed a renormalizable theory of gravity by introducing quadratic terms in curvature invariants. 
Barth and Christensen gave a detailed analysis of the one-loop divergences of fourth-order gravity theories providing 
the first general scheme of quantization of higher-order theories [1^, [5| ■ Several results followed and today it is well 
known that a renormalizable theory of gravity is obtained, at least at one-loop level, if quadratic terms in the Ricmann 
curvature tensor and its contractions are introduced [3l| . Any action, where a finite number of terms involving power 
laws of curvature tensor or its derivatives appears, is a low-energy approximation to some fundamental theory of 
gravity which, up to now, is not available. For example, String Theory or Supergravity present low-energy effective 
actions where higher-order or nonminimally coupled terms appear (4lj . 

However, if Lagrangians with higher-order terms or arbitrary derivatives in curvature invariants are considered, 
they are expected to be non-local and give rise to some characteristic length l of the order of Planck length. The 
expansion in terms of R and OR, for example, at scales larger than Iq produces infinite series which should break 
near Iq [12]. With these facts in mind, taking into account such Lagrangians, means to make further steps toward 
a complete renormalizable theory of gravity. Cosmological models coming from such effective theories have been 
extensively studied by the Noether Symmetry Approach. In particular, nonminimally coupled theories of the form 



F(<p)R+^V^V^-V(<p) 



(42) 



where F(ip) and V(tp) are respectively the coupling and the potential of the scalar field, and fourth-order theories like 

C = J—gf(R), (43) 

where f(R) is a generic function of the scalar curvature 2 . In [3 HEi Ei 51] > it was shown that asking for the existence 
of a Noether symmetry 



L X C = 0^ XL = 0, 



(44) 



where Lx is the Lie derivative with respect to the Noether vector X, it is possible to select physically interesting 
forms of the interaction potential V(ip), the gravitational coupling F(ip) and the function f(R). 

As we said above, the existence of Noether symmetries allows to select constants of motion so that the dynamics 
results simplified. Often such a dynamics is exactly solvable by a straightforward change of variables where a cyclic 
one is present. This occurrence reveals extremely useful in Quantum Cosmology since allows to formulate minisu- 
perspace models exactly solvable. Considering the above cases, i.e. scalar-tensor and fourth-order gravity, suitable 
2-dimensional configuration spaces (minisuperspace) can be achieved adopting a Friedmann- Robertson- Walker (FRW) 
metric which gives the pointlike Lagrangians. 

In the case of scalar-tensor theories, we have 



C = 6aa 2 F + 6a 2 aF - 6kaF 



V 



(45) 



2 The field equations of f(R) theories are of fourth-order in metric derivatives. They reduce to the standard second order Einstein 
equations for f(R) = R. 
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in terms of the scale factor a. 

The configuration space of such a Lagrangian is Q = {a,tp}, i.e. a two-dimensional minisuperspace. A Noether 
symmetry exists if (|44|) holds. The related Noether vector has to be 

d d d-d 
X = a—+f3—+a—+(3—, (46) 
oa dip oa dip 

where a, (3 depend on a, if. The condition (]44[) gives rise to a set of partial differential equations whose solutions are 
a, (3, F(ip) and V(ip) (see [l9[ for details). For example, a solution is 

a = -^p(s)(3 a s+1 ip m W- 1 , = p a s f m{s) , (47) 

F(ip) = D(s)ip 2 , V(ip) = Xip 2 ^ s \ (48) 

where 

(2s + 3) 2 3(5 + 1) 2.s 2 + 6.s + 3 

^ (s) -48( S +l)( S + 2)' P(S) -^7T3"' m(s) " 2 S + 3 ' (49) 

and s, A are free parameters. A suitable change of variables gives 

»^OflV Ks) , Z = -s-^prO-V 1 -^") , (50) 

PoX(s) 



where ctq is an integration constant and 



Lagrangian (143)) becomes, for fc = 0, 



where z is cyclic and 



^ = -^3- < 51 > 



£ = -/{s)w s/3 zw-\w, (52) 



2s + 3 
12t7g(s + 2)(s+ 1) 



V(*)= ^ f . ■ ow- , ^ (53) 



The conjugate momenta are 



tt z = — = 7 (s)w s/3 w = S , tt^ = — = j{s)w s/3 z , (54) 
Oz ow 



and then the WDW equation is 

[(id z )(id w )+\w 1+s / 3 ]\V>=0, (55) 

where A = j(s)X. 

The quantum version of the first of momenta (|54|) is 



-z9 2 |*>=E |*> , (56) 
so that dynamics results reduced. A straightforward integration of Eqs. ([55]) and (|56|) gives 

|* >= \n(w) > \ X (z) >oc e lS " 2 e ^- 2+s/3 , (57) 

which is an oscillating wave function. 

Analogously, in fourth-order gravity case, the pointlikc FRW Lagrangian is 

C = 6aa 2 p + 6a 2 ap - 6kap - a 3 W(p) , (58) 
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Gravitational theory 


P 


z 


w 


Scalar- Tensor 


(2s+3) 2 ,_2 
48(s+l)(s+2)^ 


/3 x(s) ^ 




Fourth order gravity s — —2 


f'(R) 


hip 


ap 


Fourth order gravity s = 


f'(R) 


a 2 


a 3 p . 



TABLE I: Summary of the changes of variables induced by the Noether symmetry for scalar-tensor and fourth-order minisuper- 
space models. The parameters as functions of s are x( s ) = — 6s/(2s+3), m(s) = (2s 2 +6s-|-3)/(2s-|-3) and l(s) = 3(s+l)/(2s+3). 



which is of the same form of (pf5|) a part the kinetic term., being p = f'(R). This is an Helmhotz-like Lagrangian |44| 
and the configuration space is now Q = {a,p}; p has the same role of the above ip. Condition (f44|) is now realized by 
the vector field 

d d d d 

X = a (a,p)—+P(a,p)— + a—+$— (59) 

The solution of this system, i.e. the existence of a Noether symmetry, gives a, (3 and W(p). It is satisfied for 

a = ot{a) , /3(a,p) = Poa s p, (60) 
where s is a parameter and /3o is an integration constant. In particular, 

s = 0^a(a) = -!ja, /3(p)=fl,j>, W(p) = W p, k = 0, (61) 



s = -2 — ► a(o) = , /3(a,p) = 4 , W(p) = W^iP 3 , V fc , (62) 
a cr 

where Wo and Wi are constants. As above, the new set of variables Q = {z, w} adapted to the foliation induced by 
X can be achieved. The procedure to achieve the solution for the WDW equation is exactly the same. The results 
are summarized in Tab.l. The solutions for the respective WDW equations are given in Tab. 2. By a rapid inspection 
of the wave functions in Tab. 2, it is clear that the oscillatory behavior, i.e. the presence of peaks, is strictly related 
to the Noether constant Eg. In other words, the Hartle criterion can be immediately implemented when conserved 
quantities are found out. 



VI. MINISUPERSPACE TOMOGRAMS IN STATIONARY PHASE APPROXIMATION 



The above results can be immediately translated in the tomographic representation. If Noether symmetries are 
present, the solutions of the WDW equation, take the following form 

m 

*(Qi,...,Qn)=5> lS ^ X (Q Z ) m<l<n (63) 



then we find, according to the results of Sec lIVl that the corresponding tomograms are 



W(Xi . . .X n ,Hx, ...,fln,Vi,... 
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Gravitational theory 


*(z,u>) 


Scalar- Tensor 




Fourth order gravity s = — 2 


e i[Siz+9feiu 2 +(3iDi/4)u) 4 ] 


Fourth order gravity s = 


el E [ z -(l/4)ln»] w l/2 z . (A]W) 


Fourth order gravity s — and w » 


e i[E z-(E /4) lnu>±Au>] 



TABLE II: Summary of the solutions of the WDW equation for the minisuperspace models. The Zo(\, w) are Bessel function. 



U=m+1 



(64) 



It is then straightforward to study the form of the tomograms satisfying the Hartle criterion. 

Let us consider first the tomograms resulting from the scalar-tensor theories. From Tab. 2, the wave function takes 
the form 



If we consider the case s = 0, the corresponding two variable tomogram is 

W(X 1 , X 2 , £tl,£t2,fl, v 2 ) oc 



Mil |/x 2 - 2A^ 2 | 



(65) 



(66) 



For the fourth order gravity models, we have considered the two cases s — — 2 and s = 0. In the first case, the wave 
function is 



and the corresponding tomogram is 



W(Xi, X 2 ,H\,H2,V\,V2) oc 



lMi| 



^ exp(?(9fc a 2 (X 2 , M2 ,^ 2 ) + 3fiq4(X 2 , // 2 , ^)) 



V9wi^2 + 18fcj/ 2 + A*2 



where ai, a 2 , 0:3 are the solutions of equation 



/i 2 



3^1^ + — + 18fc w - = 



v 2 



X 2 



v 2 



(67) 



(68) 



(69) 



ai = 



-3 • 2 1 / 3 (^ + 18fci/) 



^-81 X v 2 + 6561 X 2 ^ 4 + 2916 ^ 3 (^ + 18 fc i/) 3 ^) 



1/3 



+ 



81 X z; 2 + ^/6561 X 2 + 2916 ^ 3 (/i + 18 k v) 
3-2V3J/ 
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3 (l + i-y/3) (/Lt + 18fcz/) 

«2 = 



2 2 /3 ^_81 X ^ 2 + ^6561 X 2 ^ + 2916 ^ 3 + 18 fc ^ 



1/3 



(1 — i \/3) ^-81 Xv 2 + \j 6561 X 2 ^ + 2916 ^ (^ + 18 fc 

6-2V3!/ 

3 (l-i-v/3) (/i + 18fcz/) 



1/3 



a 3 = 



1/3 



. V3 

3 



2 2 / 3 ^-81 X v 2 + yj 6561 X 2 ^ + 2916 ^ 3 + 18 k vf^j 

(l + iV3) ^-81 Xv 2 + \j 6561 X 2 ^ 4 + 2916 ^ 3 (/i + 18 k v) 

6-2V3J, • 
The second interesting case is for s — o. When w » the wave function takes the form 

V(Z,W) = e iPo*-(E /4)lnu,±A«,] (70) 

and the corresponding tomogram is 



W(Xl,X 2 , Ml, M2, ^1,^2) oc -^-r 

M 



cxp I « 



"+ (71) 
M2) 1/2 



+ 



exp(i (-^ln/3 2 ±A/3 2 )) 



((S ^/4/3 2 2) +M2 ) 1 /2 



which can be rewritten in the form 



l/iil V((So^/4/3 2 ) + /i 2 ) 



W(X 1 ,X 2 , Ml , M2 ^ 1 ^ 2 ) oc — ( — - — T ; - ^ + (72) 



1 



+ 2 



((£ */ 2 /4/3 2 2 )+M2) 

cos(/3i - /3 2 ) 



((So^/4/3 2 ) + ^u^ 1 / 2 ((XW4/3!) + M2) 1/2 / 
where f3\ and /3 2 are the solutions of equation 

^-[-(S /4) \nw ± \w + -^w 2 - ^w] = -(So/4)- ± A + — w - — = (73) 

OW 2v 2 V 2 W V 2 V 2 



(±A - X 2 ) - ^{±\v 2 - X 2 ) 2 - SqM2^2 

A = o , 

2^2 



(±A - X 2 ) + ^J(±\is 2 - X 2 ) 2 - S M2^2 

#2 = ^ • 

2^2 
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In general, if we are able to achieve the WDW wave function in an explicit form it is always possible to construct the 
tomographic counterpart. 

At this point, we have to make a remark. If one takes the wave function ty(x) in the form of the De Broglie wave 
—=e lkx , the probability density |^(a;)| 2 is — . Also the tomogram calculated for this wave function is proportional 

y 7T Z7T 

to the value - — r. Using the tomogram for reconstructing the wave function, we get divergent integrals. Due to this 
\fi\ 

fact, one can use a "regularized" wave packet = A(x)e zkx with the amplitude A(x), e.g., of the Gaussian form 

_ 2 

~ e gx . In this case, the tomogram assumes the form of a Gaussian distribution and the density matrix for the 
pure state ip(x)ip*(x) can be reconstructed since all the integrals do diverge. The "regularized" wave packet, having 
a highly oscillatory behavior for g — > 0, corresponds to the Hartle criterion and then to the presence of Noether 
symmetries. In summary, when we have the tomograms as in the above examples containing \fi\ in the denominator, 
the reconstruction formula can be explained in the described sense. 



VII. DISCUSSION AND CONCLUSIONS 



In this paper, we discussed the realization of minisuperspace Quantum Cosmology adopting a tomographic prob- 
ability representation for the wave function of the universe. The physical meaning of the results is recovered if, in 
semiclassical WKB approximation, it is possible to select conserved quantities (stationary phases) where the wave 
function is peaked. In this case, the interpretative Hartle criterion can be applied and classical trajectories, corre- 
sponding to observable universes, are recovered. As a matter of fact the existence of Noether symmetries lead to 
oscillating components of the wave function, which present the picks required for a transition from the quantum 
cosmological states to the classical universes. In our approach, we found that a Noether symmetry So implies a factor 
of the form l/|/i| in the tomogram. Noteworthy the tomograms assume a form which is very close to classical ones, 
showing in a natural way the transition from the quantum initial stages of the universe to its classical evolution. 
This fact is crucial in our discussion since tomograms can allow, in principle, a full description of the universe from 
its initial quantum states up to the today observed. Specifically, considering the above considered models, it is easy 
to find out, either from tomograms or from Noether constants, exact classical solutions. In the case of scalar-tensor 
models, we get 

w(t) = [ht + fc 2 ] 3 /( s+3 >, (74) 
z{t) = [fcxt + fc 2 p +6 )/< s+3 ) + z , (75) 

which can be immediately translated into the original configuration space Q = {a, tp}, that is 

a(t) = a (t~t ) h ^ 7 (76) 
<p(t) = cp (t-to) q{s \ (77) 

where 

, 2s 2 + 9s + 6 , . 2s + 3 
s(s + 3) s 

Depending on the value of s, we get Friedman, power-law, or pole-like behaviors, that is all the standard classical 
cosmological behaviors. Analogously, for the fourth-order models, we have, for s = 0, 

a(t) = a e( A /6)* exp {-fe^} , (79 ) 



p(t) = p e (A/6)t exp {z x e -( 2A / 3 ) 4 } , (80) 

where ao,po and z\ are integration constants. It is clear that A plays the role of a cosmological constant and inflationary 
behavior is asymptotically recovered. For s = —2, we get power-law behaviors for a(t) and p(t). 

Such solutions, in principle, give rise to "observable universes" since the set of parameters {Hq, qg, Qm, ^a, ^k} can 
be obtained, in a standard way, from the solutions a(t), ip(t) and p(t) [45| . This fact could be extremely relevant, also 
in view of the recent observational trends which have given rise to the so called Precision Cosmology (see for example 
[H, S3j 51] ) ■ In fact, the possibility to trace back the dynamics, via tomograms, from the today observed parameters 
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up to the initial quantum conditions could be an interesting approach to formulate comprehensive cosmological models 
enclosing early and late evolution. In other words, due to this feature, the tomographic approach could be, with respect 
to other approaches included the Wigner function one, the most suitable to allow a connection between inflation and 
dark energy cosmology. 

To conclude, we point out that all the main results of Quantum Cosmology can be obtained using, arbitrarily, 
some well known formulation of Quantum Mechanics [4^ | including the density matrix and the Wigner function. 
This means that, in principle, also other approaches could be successfully used with significant results [50(. However, 
as emphasized here, the formulation based on tomographic probability distribution has some interesting properties 
which are convenient to be used as soon as the quantum-classical transitions or the comparison between classical and 
quantum pictures are relevant: in this framework, tomograms assume a fundamental role being objects which fairly 
describe states both in classical and quantum domains. 
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